Full Counting Statistics of Multiple Andreev Reflections 
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We derive the full distribution of transmitted particles through a superconducting point contact of 
arbitrary transparency under voltage bias. The charge transport is dominated by multiple Andreev 
reflections. The counting statistics is a multinomial distribution of processes, in which multiple 
charges ne (n = 1,2,3,...) are transferred through the contact. For zero temperature we obtain 
analytical expressions for the probabilities of the multiple Andreev reflections. The current, shot 
noise and high current cumulants in a variety of situations can be obtained from our result. 
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The complete understanding of the electronic trans- 
port in mesoscopic systems requires information that 
goes beyond the analysis of the current. This explains 
the great attention devoted in the last years to current 
fluctuations in these systems jlj . An important goal is to 
obtain the full current distribution. This was realized by 
Levitov and coworkers '^1, who borrowed the concept of 
full counting statistics (FCS) for photons and adapted it 
to electrons in mesoscopic systems. FCS gives the prob- 
ability P{N) that N charge carriers pass through a con- 
ductor in the measuring time. From the knowledge of 
these probabilities one can easily derived not only the 
conductance and noise, but all the cumulants of the cur- 
rent distribution. Since the introduction of FCS for elec- 
tronic systems, the theory has been sophisticated and ap- 
plied to many different contexts (see a recent review, j^])- 
In particular one of the authors and Nazarov have shown 
that, based on a Keldysh-Green's function method, one 
can calculate in a unified manner the FCS of all contacts 
involving superconducting elements 5]. 

In the context of superconductivity the basic situation, 
in which the FCS has not been yet investigated, is a point 
contact between two superconductors out of equilibrium. 
In this system the transport properties for voltages V 
below the superconducting gap A are dominated by co- 
herent multiple Andreev reflections (MAR) In these 
processes a quasiparticle undergoes a cascade of Andreev 
reflections until it reaches an empty state in one of the 
leads. Recently, the microscopic theory of MAR |^ has 
provided a new insight into this problem and has allowed 
the calculation of properties beyond the current such as 
the shot noise |^. The predictions of this theory have 
been quantitatively tested in an impressive series of ex- 
periments in atomic-size contacts [i|HjJll|. In particu- 
lar, the analysis of the shot noise |3, lul has suggested, 
that the current at subgap energies proceeds in "giant" 
shots, with an effective charge q ~ e(l -I- 2A/\eV\). How- 
ever, strictly speaking, the question of whether the charge 
in these contacts is indeed transferred in big chunks can 
only be rigorously resolved by the analysis of the FCS. 
This leads us to the central question addressed in this 



paper: what is the FCS of MAR? 

The answer, which we derive below, is that the statis- 
tics is a multinomial distribution of multiple charge 
transfers. Technically, we find that the cumulant gen- 
erating function (CGF) for a voltage V has the form 
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The CGF is related to the FCS by P{N) = 
f^_^{dx/2n) exp [S{x) — iNx]- The different terms in the 
sum in Eq. JQ) correspond to transfers of multiple charge 
quanta ne at energy E with the probability Pn{E,V), 
which can be seen by the (27r/n)-periodicity of the ac- 
companying x-dependent counting factor. This is the 
main result of our work and it proves, that the charges 
are indeed transferred in large quanta. 

To arrive at these conclusions, we consider a voltage- 
biased superconducting point contact, i.e. two super- 
conducting electrodes linked by a constriction, which is 
shorter than the coherence length and is described by a 
transmission probability T. To obtain the FCS in our 
system of interest we make use of the Keldysh-Green's 
function approach to FCS introduced Nazarov and one 
of the authors 4, 5:]. The FCS of superconducting con- 
strictions has the general form 
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Here Gi(2) denote matrix Green's functions of the left 
and the right contact. The symbol (8) implies that the 
products of the Green's functions are convolutions over 
the internal energy arguments, i.e. (Gi ^ G2){E, E') = 
J dEi Gi{E,Ei)G2{Ei,E'). The trace runs not only 
over the Keldysh-Nambu space, but also includes inte- 
gration energy. For a superconducting contact at finite 
bias voltage the CGF depends on time and Eq. ||2Jl is in- 
tegrated over a long measuring time to, much larger than 
the inverse of the Josepshon frequency. 

Let us now describe the Green's functions entering 
Eq. The counting field x incorporated into the 
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FIG. 1: Schematic representation of the MARs for BCS su- 
perconductors with gap A. We have sketched the density of 
states of both electrodes. In the upper left panel we describe 
the process in which a single electron tunnels through the sys- 
tem overcoming the gap due to a voltage eV > 2 A. The other 
panels show MARs of order n — 2,3, 4. In these processes an 
incoming electron at energy E undergoes at least n — 1 An- 
dreev reflections to finally reach an empty state at energy 
E + neV. In these MARs a charge ne is transferred with a 
probability, which for low transparencies goes as T". At zero 
temperature they have a threshold voltage eV — 2A/n. The 
arrows pointing to the left in the energy trajectories indicate 
that a quasiparticle can be normal reflected. The lines at en- 
ergies below E and above E + neV indicate that after a detour 
a quasiparticle can be backscattered to finally contribute to 
the MAR of order n. 



matrix Green's function of the left electrode as follows 
Gi (x, t, t') = e-^x^^/^Gi {t, t')e*x^^/2 (3) 

Here Gi{t,t') is the reservoir Green's function in the ab- 
sence of the counting field and tk = 173 fa a matrix in 
Keldysh(")-Nambu(~) space. We set the chemical po- 
tential of the right electrode to zero and represent the 



Green's functions by Gi{t, t') = e 



and G2{t, t') = Gs{t—t'). Here, we have not included the 
dc part of the phase, since it can be shown that it drops 
from the expression of the dc FCS at finite bias. Gg is 
the Green's function of a superconducting reservoir (we 
consider the case of a symmetric junction), which reads 
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tions of the leads and f{E) is the Fermi function. Ad- 
vanced and retarded functions in have the Nambu- 
structure R{A) = g^'^f^ f^'^fi fulfilling the normal- 
ization condition + = 1. They depend on energy 
and the superconducting order parameter A. 

In Eq. (0) the matrix appearing inside the logarithm 
has an infinite dimension in energy space. In the case 
of N-N or N-S contacts such a matrix is diagonal in this 
space, which makes almost trivial the evaluation of the 
FCS. In the S-S case at finite bias this is no longer true, 
which introduces an enormous complication. 

We now tackle the problem of how the functional con- 
volution in Eq. ||2J) can be treated. The time-dependence 
of the Green's functions leads to a representation of the 
form G{E,E') = Y.nGo,n{E)S{E - E' + neV), where 
n = 0,±2. Restricting the fundamental energy inter- 
val to E — E' E [0, eV] allows to represent the con- 
volution as matrix product, i.e. (Gi (E) G2){E,E') 

{GiG2)n.m{E,E') — J2ki^l)n^k{E,E'){G2)k.m{E, E'). 

Writing the CGF as Six) = {to/h)TTlnQ, where Q = 
1 + {VT/2) (Gi(x) - G2) El- The trace in this new 

representation is written as /g'^^ di^X^n ^i' In (Q)„„- In 
this way the functional convolution is reduced to matrix 
algebra for the infinite-dimensional matrix Q. Still, the 
task to compute TrlnQ is nontrivial. However, noting 
that TrlnQ = IndetQ, it is obvious at this stage that 
det Q has the form of a Fourier series in x, which allows 
us to write the CGF as follows 
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Keeping in mind the normalization S'(O) = 0, it is clear 
that one can rewrite this expression in the form antic- 
ipated in Eq. (1), where the probabilities are given by 
PniE, V) = P:,iE, V)l En=-oo KiE, V). Of course, one 
has still to extract the expression of these probabilities 
from the determinant of Q, which is a non-trivial task. 
It turns out that Q has a block-tridiagonal form, which 
allows to use a standard recursion technique. We define 
the following 4x4 matrices 
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Here i?(^)(£') are retarded and advanced Green's func- 



wheTeQn,miE) = QiE+neV, E+meV). With these defi- 
nitions, det Q is simply given by det Q = Ilj^-oo "^^^ ^'^r 
In practice, det F„ = 1 if \n\ :§> A/|eV"|. This reduces the 
problem to the calculation of the determinants of 4 x 4 
matrices. 

In the zero-temperature limit one can work out this 
idea analytically to obtain the following expressions for 
the probabilities 
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Here, we have used the shorthand g^'^'{E) = g^'^{E + neV), and defined 
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where a = i?, A, Kn,m — (Iljli <i"3t ^n-2j)(njli "^^t Fm+2j) and the diS^erent functions can be expressed as foUows 
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Notice that, since at zero temperature the charge only 
flows in one direction, only the Pn with n > sur- 
vive. It is worth stressing that the full information about 
the transport properties of superconducting point con- 
tacts is encoded in these probabilities. Let us remark 
that Pn{E,V) are positive numbers bounded between 
and 1. Although at a first glance they look com- 
plicated, they can be easily computed and provide the 
most efficient way to calculate the transport properties 
of these contacts. In practice, to determine the functions 
B:^-^ and detK, one can use the boundary condition 
B:^'^ = detK = 1 for \n\ > A/\eV\. For perfect trans- 
parency (T — 1) the previous expressions greatly simplify 
and the probabilities Pn{E, V) can be written as 
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FIG. 2: Current contribution of processes n = 1, 10, from 
right to left, as a function of voltage for BCS superconductors 
of gap A. The different panels correspond to different trans- 
missions. Notice the logarithmic scale in the upper panels. 



where a{E) is the Andreev reflection coefficient defined 
as a{E) = -if'^iE)/ [l + g'^iE)] , and a„ = a{E + neV). 

In view of Eqs. H7I9|I the probabilities P„ can be in- 
terpreted in the following way. P„ is the probability of 
a MAR of order n, where a quasiparticle in an occupied 
state at energy E is transmitted to an empty state at 
energy E + neV. The typical structure of the leading 
contribution to this probability consists of the product 
of three terms. First, Jq gives the probability to in- 
ject the incoming quasiparticle at energy E. The term 
rifc^i (^/4)|/fc^P describes the cascade of n — 1 Andreev 
reflections, in which an electron is reflected as a hole and 
vice versa, gaining an energy eV in each reflection. Fi- 
nally, J„ gives the probability to inject a quasiparticle in 
an empty state at energy E + neV. In the tunnel regime 



Pn{E,V) = (rV4"-i)poP„nfc=i'l/fe^P, PiE) being the 
reservoir density of states. This interpretation is illus- 
trated in Fig. ^ where we show the flrst four processes 
for BCS superconductors. The product of the determi- 
nants in the expression of (see Eq. (fTj)) describes the 
possibility that a quasiparticle be reflected and make an 
excursion to energies below E or above E + neV [l3| . 
In the tunnel regime this possibility is very unlikely and 
at perfect transparency is forbidden. As can be seen in 
Eq. UlUI) . for T — 1 the quasiparticle can only move up- 
wards in energy due to the absence of normal reflection. 

From the knowledge of the FCS one can get a deep 
insight into the different transport properties by ana- 
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FIG. 3: Second and third cumulant at zero temperature for a 
contact between BCS superconductors. Both are normahzed 
to the first cumulant (the average current). The transmissions 
are indicated in the plots. 



lyzing the role played by every process. For instance, 
in Fig. 121 we show the contribution to the dc current of 
the individual processes, i.e. /„ = {2e/h) J dE nP„, for 
the case of BCS superconductors of gap A. In this case 



iA/ [{E =F iSf - A^] , where S = 0+, and g 



A,R 



follows from normalization. As can be seen in Fig. El a 
MAR of order n has a threshold voltage eV — 2A/n, 
below which it cannot occur. The opening of MARs at 
these threshold voltages is the origin of the pronounced 
subgap structure visible in the different transport prop- 
erties (see Fig. O. Notice also that at low transmission 
the MAR of order n dominates the transport for voltages 
[2A/n, 2A/(n — 1)], while at high transparencies several 
MARs give a significant contribution at a given voltage. 
This naturally explains why the effective charge is only 
quantized in the tunnel regime 0, Q| . 

From the CGF one can easily calculate the cumulants 
of the distribution and in turn many transport properties. 
Of special interest are the first three cumulants Ci = N, 
C2 = {N -N)'^ and C3 = {N -7^)^, which correspond 
to the average, width and skewness of the distribution, 
respectively. From the fact that the FCS is a multinomial 
distribution, it follows that at zero temperature these 
cumulants can be expressed in term of the probabilities 
Pn{E, V) as Cn{V) = dE CniE, V), where 

00 00 

Ci{E, v) = J2 ; C2{E, v)^J2 - ^1 ^) 

n— 1 n— 1 

00 

CaiE, V)^Y1 - [^1 ^) + 3^2(^^, V)] . 

n=l 

The first two cumulants are simply related to the dc 
current, / — {2e/h)Ci, and to the zero-frequency noise 
Si — [Ae^ /h)C2- In Fig. we show C2 normalized by 
Ci, which reproduces the results for the shot noise re- 
ported in the literature . In this figure we also show C3. 



This cumulant determines the shape of the distribution, 
and it is attracting considerable attention El [la] be- 
cause it contains information on nonequilibrium physics 
even at temperatures larger than the voltage. As seen in 
Fig. O at low transmissions C3 = g^Ci, where q{V) — 
1 + Int(2A/ey) is the charge transferred in the MAR 
which dominates the transport at a given voltage. This 
relation is a striking example of the general relation con- 
jectured in Ref. and it is simply due to the fact 
that the multinomial distribution becomes Poissonian in 
this limit. For higher transmissions this cumulant is 
negative at high voltage as in the normal state, where 
C3 = T(l - r)(l - 2T), but it becomes positive at low 
bias. This sign change is due to the reduction of the 
MAR probabilities at low voltage. After the sign change 
there is a huge increase of the ratio C3 /Ci , which is a sig- 
nature of the charge transfer in large quanta. Finally, at 
T = 1 the cumulants (C„ with n > 1) do not completely 
vanish due to the fact that at a given voltage different 
MARs give a significant contribution, and therefore their 
probability is smaller than one (see Fig.|2td)). 

In summary, we have demonstrated that in supercon- 
ducting contacts at finite voltage the charge transport 
is described by a multinomial distribution of multiple 
charge transfers. This proves that in the MAR processes 
the charge is indeed transmitted in large quanta. We have 
obtained analytically the MAR probabilities at zero tem- 
perature, from which all the transport properties are eas- 
ily computed. Our result constitutes the culmination of 
the recent progress in the understanding of MARs, which 
are a key concept in mesoscopic superconductivity. 
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